Phase separation in the vicinity of "quantum critical" doping concentration: 
implications for high temperature superconductors 
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A general quantitative measure of the tendency towards phase separation is introduced for systems 
exhibiting phase transitions or crossovers controlled by charge carrier concentration. This measure 
is devised for the situations when the quantitative knowledge of various contributions to free en- 
ergy is incomplete, and is applied to evaluate the chances of electronic phase separation associated 
with the onset of antiferromagnetic correlations in high-temperature cuprate superconductors. The 
experimental phenomenology of lanthanum- and yittrium-based cuprates was used as input to this 
analysis. It is also pointed out that Coulomb repulsion between charge carriers separated by the 
distances of 1-3 lattice periods strengthens the tendency towards phase separation by accelerating 
the decay of antiferromagnetic correlations with doping. Overall, the present analysis indicates that 
cuprates are realistically close to the threshold of phase separation — nanoscale limited or even 
macroscopic with charge density varying between adjacent crystal planes. 

PACS numbers: 



I. INTRODUCTION 



Phase transitions or sharp crossovers induced by con- 
tinuous change in chemical composition have natural ten- 
dency to cause phase separation. Phase separation oc- 
curs, when the free energy of the most stable homoge- 
neous state has concave dependence on the concentration 
of particles (i.e. this dependence is curved or peaked up- 
ward). The concave behavior is particularly likely in the 
course of a phase transition or a crossover, because the 
free energy switches between two different dependences 
characteristic of adjacent phases. The steeper is this 
switching, the likelier is the onset of concave behavior. 
In the case of electronic subsystems in solids, a strong 
Coulomb repulsion raises the threshold of phase separa- 
tion and also, if the homogeneous electronic state be- 
comes unstable, limits the size of single-phase regions to 
the nanometer scale along at least one of the spatial di- 
rections. In the present work, we introduce a quantitative 
measure of the above tendency towards phase separation 
in the vicinity of phase transitions and crossovers for the 
situations when the knowledge of various contributions 
to free energy is incomplete. We justify this measure in 
a very general context, and then use it to analyze the 
chances of electronic phase separation in high tempera- 
ture cuprate superconductors. 

Macroscopic phase separation in cuprates had been 
observed in LaCuC^+a T] , where it occurs because of 
the high mobility of intercalated oxygen that maintains 
charge neutrality within each phase. In other cuprate 



families, where dopant ions are immobile, the experimen- 
tal situation is complicated by the possibility that phase 
separation may be limited by Coulomb interaction to the 
difficult-to-access scale of a few nanometers, where in- 
homogeneities can also fluctuate in time. The direct ex- 
perimental evidence for static nanoscale modulations has 
been reported so far only for a rather limited subset of 
cuprates 0, S,i,[E[E 0,[|. 

Theoretically, the general tendency towards phase sep- 
aration in cuprates and cuprate-related models is exten- 
sively discussed in the literature 
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but any conclusion about the actual presence or absence 
of phase separation suffers from uncertainties associated 
with the gross simplifications in the models. The impor- 
tant advantage of the phase separation measure intro- 
duced in this work is that it is directly applicable to real 
materials. At the same time, our analysis also provides a 
useful perspective to the results of various model studies. 

Cuprates can exhibit phase separation for two related 
reasons: the onset of antiferromagnetic (AF) correla- 
tions and the onset of Mott insulating gap. In this 
work, we focus on the former and arrive at the conclu- 
sion that cuprates are very close to the phase separa- 
tion threshold and can realistically phase separate. The 
The experimental results on the doping dependence of 
the overall intensity of antiferromagnetic fluctuations in 
YBa 2 Cu306+a;fll| and La9--rSr. r CuO,i[32| constitute im- 
portant input into our analysis. We also propose how the 
phenomenology of phase separation in LaCu04+,s[l[ can 
be used to obtain information about the properties of 
"generic" cuprates. Finally we point out that Coulomb 
repulsion between charge carriers separated by the dis- 
tances of 1-3 lattice periods strengthens the tendency to- 
wards phase separation by accelerating the decay of AF 
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with two constraints 



h + h = 1, 



(3) 
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FIG. 1: (Color online) Temperature- vs. -concentration phase 
diagram for homogeneous phases. The solid line can represent 
a phase transition of any order or a crossover. 



correlations with doping. 



II. PHASE TRANSITIONS AND PHASE 
SEPARATION 

A. Free energy 

We analyze the dependence of free energy of energet- 
ically most stable homogeneous state as a function of 
charge carrier concentration. Negative second derivative 
of this dependence (equivalent to negative curvature or 
negative inverse compressibility) will be considered as a 
sign of phase separation. Motivated by the physics of 
cuprates, we consider the "homogeneous" temperature- 
vs. -concentration phase diagram looking as shown in 

Fig. ni 

The total free energy of the system i<t t as a function 
of dimcnsionless charge carrier concentration x can be 
decomposed as shown in Fig. [2l namely: 



F tot {x)=F Q (x)+F n (x), 



(1) 



where F (x) is the free energy of a non-ordered state 
(e.g. non-magnetic Fermi liquid), and F^x) is the part 
of free energy associated with the onset of a more or- 
dered state below the transition or crossover at x = x c0 
(e.g. AF Neel phase or some kind of electronic liquid 
with strong AF correlations like resonating valence bond 
(RVB) state jH). 

Let us now introduce the Maxwell construction by as- 
suming that we are dealing with an idealized neutral sys- 
tem, where particles separate into two macroscopic re- 
gions: fraction j\ of them forms regions with the concen- 
tration x\, while fraction f 2 goes into the concentration 
X2 ■ The free energy -Fsep of this phase separated mix is 
given by equation 



F f 



sep 



fiF t ot(xi) + f2F tot (x 2 ) 



(2) 



(4) 



where x a v is the average concentration set by the doping 
level. The minimization of -Fsep with respect to /i, f 2 , 
x\ and X2 under conditions (j3j [4]) gives: 



e^tot 



dx 



dF tot 



dx 



-Plot (£2) ~ ^totOi) 



X2 - Xi 



(5) 



The derivatives ^ot are ^ f coursej the chemical poten- 
tials of the two phases. 

The above condition has transparent geometrical in- 
terpretation — Maxwell construction — illustrated in 
Fig-[2Ib): The line connecting points (xi, -Ftot( x i)) an d 
(x2, -Ftot^)) should be tangential to the curve i 7 tot( a; ) 
at the both points. The energy of the phase separated 
state is then a point on this line at x = x av . The nec- 
essary condition for Eq.((5]) to correspond to a minimum 
rather than a maximum is to have a region of negative 
curvature somewhere between X\ and X2. 

If i*tot ( x ) an d its first derivative are continuous every- 
where, then there exist three different regions of carrier 
concentrations: 

1) Miscibility gap x m < x < x c o, where the curvature 
is negative, and hence, the uniform composition locally 
unstable towards spinodal decomposition. 

2) Metastable regions x\ < x < x m and x c q < x < X2, 
where the homogeneous state is locally stable due to the 
positiveness of the curvature, but globally unstable, be- 
cause the fully separated inhomogeneous state has lower 
energy. 

3) Stable regions x < xi and X2 < x. 

If i*i(x) exhibits a cusp 34], then one of the above re- 
gions can shrink to a point. 

Now we return to the discussion of real systems de- 
scribed by free energy (JT|). We assume that the non- 
ordered state does not phase separate on its own, and, 
therefore, the curvature of Fq(x) 7 which we define as 



K 



2 dx 2 



(6) 



is positive everywhere. [This and other curvatures intro- 
duced below are implied to be x-dependent.] The onset 
of the more ordered state should necessarily lower the 
energy of the system, which, as obvious from Fig. [2fa), 
implies that F v (x) should have either intervals of finite 
negative curvature, or a point of infinite negative curva- 
ture (upward cusp). We define the "negative" curvature 
of F J1 (x) as 



1 d 2 F n 
2~dx T 



(7) 
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FIG. 2: (Color online) Sketches of free energy as a function of 
charge carrier concentration for a phase transition accompa- 
nied by phase separation, (a) Free energy of non-ordered state 
Fo(x), and free energy gain F v (x) due to the phase transition 
at Xco- (b) Total free energy F^ Q ^(x) (thick line) equal to 
the sum of Fo(x) and F v (x) sketched in (a), and the Maxwell 
construction. 



The curvature of the total homogeneous free energy can 
now be expressed as 



tot = 77 



1 d 2 F tot 

2 dx 2 



= Kn — K„ 



(8) 



A negative value of -Ktot is necessary but not sufficient 
condition for the onset of phase separation. Once an ex- 
cess charge carrier concentration Ax starts building in 
some part of the sample, there is an additional energy 
cost AF(Ax) associated with the Coulomb repulsion of 
uncompensated charge, strain energy and gradient en- 
ergy. When dopant atoms are immobile, we expect the 
main additional contribution to come from Coulomb en- 



ergy and denote the positive curvature associated with 
this contribution as 



K, 



Coul 



1 d 2 AF 

2 dAx 2 ' 



(9) 



The condition for the onset of phase separation then be- 
comes: 



K V >K + K t 



Coub 



(10) 



and the Maxwell construction bounds on the miscibility 
gap and the metastable regions should also be modified 
accordingly. 

In subsection lll Cl we quantify how large K v can be for 
the situations when only the values of -^(0) and x c q are 

dF (x) 

known, with or without additional knowledge of — 4^ 
at x = x c o and/or x = 0. 



B. Different scenarios of phase separation 

The fact that the system phase separates is sometimes 
perceived as implying first order phase transition. How- 
ever, in the case of volume constrained systems, e.g. elec- 
trons on the lattice, the relation between phase sepa- 
ration and the order of phase transition becomes less 
straightforward. In contrast to the more familiar con- 
straint of constant pressure, which allows the phase coex- 
istence only at a fixed temperature of a first order phase 
transition, the constraint of constant volume allows the 
coexistence of two phases in a range of temperatures (35|. 
Such a property, in turn, defeats the classification of the 
resulting behavior in terms of the jumps of thermody- 
namic derivatives, because, when the temperature of the 
system decreases past the phase separation threshold, the 
volume fraction of the second phase increases gradually, 
which normally implies no jump in entropy and no la- 
tent heat, but leads to a jump in the specific heat like in 
a second order phase transition. Yet, at the phase sep- 
aration temperature, the homogeneous state usually re- 
mains metastable, which means, that, in order to become 
phase separated, the system has to overcome a nucleation 
threshold. 

In our classification, we distinguish between the types 
of "would-be" homogeneous phase transitions, which are 
behind the onset of phase separation. Eventual phase 
separation is compatible with the possibilities that the 
line in Fig. [T] represents a phase transition of first order 
(Fig. [3|) — always exhibiting a cusp, second or higher or- 
ders (Fig. [3J) — with or without cusp dependently on the 
critical exponents, or it may also represent a crossover 
(Fig. 2]) . These are distinct scenarios of phase separation 
with distinct properties. Before proceeding further, the 
readers may choose to read the discussion of these scenar- 
ios given in Appendix [Al The part of this Appendix on 
the second order phase transitions exemplifies the gen- 
eral constraint on the negative curvature to be derived 
in the next subsection. 
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FIG. 3: (Color online) Sketch of a first order phase transi- 
tion. The system switches from a potential minimum Fa(x) 
associated with state A to the potential minimum Fb{x) as- 
sociated with state B. This transition necessarily leads to an 
upward cusp. 



energy 




FIG. 4: (Color online) Sketch of free energy gain F v (x) asso- 
ciated with a crossover (solid line) . Dashed line illustrates the 
graphical meaning of convent ion (|1 3 [I for fixing the location of 
the crossover point x c q. 



C. Rigorous constraints on the negative curvature 

In this subsection, we consider the following problem 
and its generalizations. 

Problem A — Phase transition: Let us assume, that we 
can isolate certain contribution F v (x) to the free energy 
of the system as a function of charge carrier concentra- 
tion x. About this contribution, we only know that, at 
x = 0, F v (0) = —F n o < and, for carrier concentrations 
greater than some given positive value x c0 , F 71 (x) = 0. 
[The sketch for F v (x) is given in Fig. H^a).] The latter 
condition implies that F v (x c o) = and F'(x c q + 0) = 0. 
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FIG. 5: (Color online) Parabola (solid line) minimizing the 
value of the maximum negative curvature for the boundary 
conditions of Problem A, and some obvious alternatives to 
this parabola (dashed lines). The straight-line alternative 
ends in an upward cusp at x = x c o, which means a cusp 
with infinite negative curvature. Two other alternatives have 
intervals with finite negative curvatures larger than that of 
parabola. 



(Here and everywhere, F'(x) refers to the first derivative 
with respect to x, and F"(x) refers to the second deriva- 
tive.) Any path leading from F v (x c o) to F n (0) should 
have a region of negative curvature (or a single point 
with infinite negative curvature). Let us further assume 
that, for every possible path F r/ (x) satisfying the above 
conditions, we can find the value of the maximum neg- 
ative curvature (max[-if^'(a;)]) in the interval [0,x c o]. 
The question we ask is: What is the minimum possible 
value of max[— iF^'(x)]? We denote this value as K v m- 
A rigorous solution of this problem is presented in Ap- 
pendix [B] The result is 



F, 



77O 



(11) 



which is equivalent to the expression (|A9[) obtained in 
Appendix [X] in the framework of Landau expansion for 
a second order phase transition. Curvature (|lip corre- 
sponds to the parabola shown in FigJ5] Other sketches 
in the same figure represent obvious alternatives to the 
parabolic dependence. These alternatives exhibit cither 
regions of larger negative curvature, or a point of in- 
finitely negative curvature (cusp) at x = x c o — both 
are the possibilities that are stronger in favor of eventual 
phase separation. 

Constraint (fTTj) will play the central role in our subse- 
quent analysis of phase transitions and crossovers. 

Problem A can be straightforwardly generalized to the 
following one. 
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Problem Al Crossover: Assume that 

F v (x c0 ) = -F nC and F^(x c0 ) = F^ c < (F v0 - F vC )/x c0 , 
where F v c and F' c are two new constants. All other 
conditions are the same as in Problem A. 

The solution of Problem Al is (see Appendix [BJ 



K 



F v o — F n c — F„ 



(12) 



In the case F' c > (F v0 ~ F n c)/XcQ, there are paths be- 
tween x = and x — x c q, which do not exhibit negative 
curvature. 

Expression (fT2")) is useful, when one needs to estimate 
the error associated with applying the basic expression 
(|llj) to a crossover, in which case the choice of x c o may 
be somewhat arbitrary. That choice, however, can be 
made unambiguous by requiring, e.g., that 



F, 



F^ c x c0 = O.l^o. 



(13) 



In this case, the error due to the application of expression 
(fTT|) is limited to 10 per cent. 

In some physical situations, the value of F'(0) may 
also be known. This knowledge can then be exploited as 
follows. 

Problem A2: F' (0) = F' , where F' is a new constant. 



All other conditions are the same as in Problem Al. 

The solution of Problem A2 is given in the AppendixlBl 
Here we only give the result: 

1. If F^ c > (F v o-F n c)/x c o > F' nQ , then the paths be- 
tween x = and x — x c q can have no negative curvature 
regions or points. 

Otherwise: 

2a. If 



2(F v0 - F nC ) 



XcO 



(14) 



Then the minimum value of the maximum negative cur- 
vature is still given by Eq. (fl2)) . 

2b. If condition (flJl is not fulfilled, then 



K, 



F 



F v o + F^qXcO 



(15) 



1 c0 



be large (in particular, in hole-doped cuprates), because 
the disappearance of the static order is due to the the 
loss of three-dimensional AF correlations between dif- 
ferent Cu02 planes. In the paramagnetic state, two- 
dimensional AF correlations should still carry large mag- 
netic energy in a broad parameter region near the transi- 
tion line. The subsequent disappearance of magnetic en- 
ergy associated with AF correlations can proceed either 
as a smooth crossover, or through one or several phase 
transitions involving possibly spin liquid states such as 
RVB. 

In order to impose the negative curvature constraint 
(llip on the doping evolution of magnetic energy, we do 
not need to know exactly what happens as the doping 
level increases. All we need to know is the easily acces- 
sible value of exchange energy at zero doping (F v q) and 
the value of critical doping x c q, where condition (jlip is 
applicable. 



B. Decomposition of free energy 

In Eq.([T]), we assign Fq(x) to be the lowest possible 
energy of a non-magnetic state at a given charge carrier 
concentration, and F v (x) to be the of the energy associ- 
ated with the onset of AF correlations both static and 
dynamic. The latter includes both the AF exchange en- 
ergy as such plus the change of kinetic, Coulomb and all 
other contributions to the total energy caused by AF cor- 
relations. With this definition, once the AF correlations 
disappear, the non-magnetic state of the lowest energy 
is the actual physical state of the system. At the dop- 
ing concentrations corresponding to non-zero AF correla- 
tions, the notion of the lowest energy non-magnetic state 
should be viewed as a variational ansatz, where the AF 
correlations are completely suppressed at and below typ- 
ical exchange frequencies. It should be analogous to the 
paramagnetic state at temperatures much higher than 
exchange coupling. At half doping (x = 0), the optimal 
non-magnetic state should be a Mott insulator, where 
localized spins have completely random orientations. 



Estimates of positive and negative curvature 
contributions to free energy 



III. ANTIFERROMAGNETIC CROSSOVER IN 
CUPRATES 

A. Preliminary remarks 

In this part, we focus on the doping dependence of 
the energy of AF correlations in the cuprates, where the 
exchange coupling between Cu spins leads to establish- 
ing AF order at half- filling (x = 0). As doping level x 
increases, static AF order disappears through a second 
order phase transition. However, the magnetic energy 
change associated with this phase transition should not 



In this subsection we estimate curvatures Kq, i^Coul 
and K v for the paramagnetic, Coulomb and antiferro- 
magnetic contributions to free energy, respectively. We 
are looking for conditions sufficient for phase separation 
as follows from inequality (jTTJjl . Therefore, at every step 
we try to overestimate the factors opposing phase sep- 
aration, and underestimate factors leading to it. [The 
constraint (|11[) works this way] 

Everywhere below, the free energy and the correspond- 
ing curvature are evaluated per in-plane Cu atom. Occa- 
sionally, we interchange the terms "doping" and "charge 
carrier concentration" . 
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1. Paramagnetic curvature 

The positive curvature of Fq(x) can be estimated as 
that of 2D Fermi-liquid 



^ 2 (l + /o) 



2m*a 2 



(16) 



where m* is the effective mass and f§ the Landau 
Fermi-liquid parameter (36|. Due to two-dimensionality, 
expression (|16p is helpfully independent of the Fermi- 
momentum and, hence, of the charge carrier concentra- 
tion. Therefore, the above estimate would apply both to 
the conventional case, when all charge carriers form the 
Fermi sea, and, to an alternative case, when only doped 
charge carriers form a Fermi-liquid in the [remnants of] 
lower Mott-Hubbard band. 

The effective mass can be estimated as 



3ft 2 7 



(17) 



where fcs is Boltzmann coefficient, and 7 is electronic 
specific heat coefficient per in-plane Cu atom for over- 
doped samples at temperatures high enough to de- 
stroy magnetic correlations. (Here, we ignore the low 
temperature pseudogap behavior in underdoped sam- 
ples as the pseudogap itself can be the consequence of 
nanoscale-limited phase separation.) For the estimate, 
we use 7 = 1 mJ/gat.K 2 [converted per in-plane Cu 
atom as required by Eq. (fT"7|) ]. This value is consistent 
with experiments on L^-zSr^CuO^, YBa2Cu306+£c and 
Bi2Sr2CaCu208+a;[33|- It leads to m* w 4m e , where m e 
is the mass of free electron. We are not aware of ex- 
periments, which would allow one to access the value of 
/q directly. It is, however, reasonable to assume that 
/o ~ fi ~ ^- H ere j /1 is another Fermi-liquid parameter 
determined by the effective mass according to relation 
m* = (1 + \fl)mi, w (1 + i/f)m e , where mj, is the 
band mass. [In normal He-3 at zero applied pressure, 
/ S = 10.07 and f( = 6.04^.} As long as f[ remains 
sufficiently large, assumption /q ~ f( guarantees that 
the resulting value of Kq equals approximately twice the 
value for the free electron gas . Our final estimate in sub- 
section IIII C 41 will include allowance for a factor-of-two 
uncertainty in the above value. 

The right-hand side of Eci. (fT6|) can be rewritten as 
(1 + /q)/(2!/), where v is the density of electronic states 
around the chemical potential. As such, this formula 
has a broader range of applicability: v can be extracted 
from the specific heat measurements as 37/(7rfcs) 2 irre- 
spectively of whether the substance measured is a Fermi- 
liquid, and parameter /q characterizes the interaction be- 
tween a newly added particle and the particles already 
present in the system. 

The lowest energy non-magnetic state may or may not 
be a Fermi-liquid, especially around zero doping, where 
the onset of Mott insulating behavior (charge-transfer 
gap), the presence of random paramagnetic background, 



and the localization effects due to the random Coulomb 
potential of dopant atoms can drastically modify the be- 
havior of charge carriers. The resulting metal-insulator 
transition is a fascinating subject on its own (see e.g. 
[T9L [39l| ^ . However, we feel, that, this transition (pre- 
sumably around x w 0.05) can only further contribute 
to the tendency toward phase separation and, thereby, 
further strengthen the conclusions of this work. 



2. Coulomb curvature in the "lasagna" scenario 

When the phase separation of 2D electronic systems 
is analyzed, one frequently mentioned picture is that 
of Coulomb- frustrated phase separation [3, EH, [H, HsJ. 
Many researchers expect that Coulomb-frustrated in- 
plane phase separation can explain the formation of 
stripes, checkerboards or other nanoscale-limited inho- 
mogeneous structures possibly existing in cuprates. 

In the present work, however, we explore a different 
possibility, which we call the "lasagna scenario". Ac- 
cording to this scenario, the three dimensional layered 
system phase separates into two-dimensional macroscopic 
regions of positive and negative charge lying on the top of 
each other in adjacent layers and having equal in absolute 
value and opposite in sign densities of uncompensated 
charge. These regions effectively screen each other over 
the distances larger than the separation between layers. 
This kind of screening does not require a perfectly peri- 
odic charge order along the c-axis like (H 1 1 1 — ), 

where + or — represents the sign of the charge. (Pe- 
riodic arrangement would imply an observable doubling 
of the c-axis period.) A more disordered arrangement, 
where each positive layer has at least one adjacent neg- 
ative layer and vice versa [e.g. (H 1 h H 1 — )] 

would be sufficient for the present scenario. 

The Coulomb cost of this phase separation scenario can 
be estimated as that of Coulomb interaction between un- 
compensated charges within radius equal to the distance 
between the CuC>2 planes: 



K Cou\ 



2ire 2 R s 



(18) 



where e is the charge of electron, R s is the distance 
between Cu02 planes, and e is the dielectric constant 
estimated as 30. Here we sidestep the discussion of 
multi-layered compounds and consider only single-layer 
cuprates. 

The lasagna scenario is worth exploring seriously, since 
it may well be realized in cuprates, and this would in- 
fluence strongly the interpretation of many experiments. 
At the same time, this scenario is easier to treat theoret- 
ically, since one does not need to deal with the gradient 
terms and the strongly fluctuational nature of nanoscale 
inhomogeneities. The Coulomb curvature (fl"8|) for our 
scenario is comparable with and can easily be lower than 
that of nanoscale limited phase separation. Absence of 
the nanoscale gradients is another energetic advantage 
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of this scenario. Where it may loose to the nanoscale- 
limited case is in the energy gained from phase separa- 
tion. As follows from Maxwell construction [Fig. EJb)], 
the charge densities of the oppositely charged regions in 
the phase separated state are, in general, not equal to 
each other. Imposing the equal charge condition means 
that the energy gain due to phase separation is less than 
optimal. Another extra cost of the present scenario is 
that of the suppressed interlayer hopping, which, how- 
ever, should be small. 

If the lasagna scenario can be proven to be energeti- 
cally advantageous over the homogeneous state, then it 
is guaranteed that the system is unstable towards phase 
separation — either according to the lasagna scenario, or 
with in-plane nanoscale inhomogeneities, if they lower the 
energy even further or better accessible kinetically. The 
nanoscale-limited phase separation can also take place as 
a second stage, after charge density becomes different in 
adjacent planes according to the lasagna scenario. 
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3. Antiferromagnetic curvature 

Here we apply constraint (fTTjl . We estimate the value 
of F v o in the approximation of staggered spin polariza- 
tions as 

F v o = \ JN nn s 2 = i J, (19) 

where J is the exchange coupling constant, N nn = 4 is 
the number of the nearest neighbors on the square lattice, 
and s = 1/2 is the value of Cu spin. 

The constraint (jlip with definition (|13[) for x c q now 
implies that in the interval < x < x c o there exist values 
of x such that 

K v {x) > K v m(x c o) = -4r, (2°) 

subject to a 10 percent uncertainty. The persistence of 
magnetic correlations with doping is sometimes viewed 
as indicating their strength and therefore favoring phase 
separation. The constraint (|2T)]) . however, suggests oth- 
erwise. 

The identification of x c q is the subject of the rest of 
this paper. This parameter is difficult to pinpoint, but 
the focus on finding it constitutes a new and, presumably, 
productive approach to the subject of phase separation 
in cuprates. The effectiveness of this approach is related 
to the fact that K n u depends very steeply on x c o- The 
readers should appreciate that K v m for x c q = 0.14 is 4 
times larger than for x c q = 0.28 and 4 times smaller than 
for x c o = 0.07. Therefore, as long as other contributions 
to the curvature are known within factor of two, and x c q 
is known with the accuracy of ±0.05, one can form a good 
judgment about the chances of phase separation. 

Both experiments and numerical studies seem to agree 
that there are no noticeable AF correlations above x = 
0.5. Further constraining x c q requires a more detailed 



FIG. 6: (Color online) Estimates for negative (K v ) and pos- 
itive (Ko + ^"coul) contributions to energy curvature per one 
in-plane Cu. Negative curvature is plotted as a function of 
unknown critical concentration x c o- Solid lines represent the 
estimates by formulas ((20}, (|16[) and (|18[1 with the numbers 
given in the text. Shaded areas around the lines cover the re- 
gions of the factor-of-two uncertainty for the above estimates. 



analysis of the experimental phenomenology and theo- 
retical scenarios. This will be done in Sections IIII Dl and 
MI El respectively with the end result that both of these 
analyses point to x c q w 0.14. 



4- Combining the estimates 

Before trying to pin down the value of x c q, it is useful 
to plot the bound on the value of as a function of x c q 
against the bounds on Kq + Kq ou [ thus comparing the 
two sides of condition (flU)) for phase separation. Such a 
plot is shown in Fig. [Sj 

The horizontal line corresponds to the value of Kq + 
-^Coul estimated according to formulas (flUj) and (fig]) 
with the following values of parameters: a = 4 A, /q = 6, 
to* = 4to c , R s = 6 A, e = 30. According to this esti- 
mate K = 1.35 eV and Kq ou i = 0.57 eV [both per 
in-plane Cu] . The shadowed stripe around the above line 
covers the region of the factor-of-two uncertainty in the 
estimated value. 

The line corresponding to K ri (x c o) is obtained using 
K v m(x c o) given by Eq.® with J = 125 meV. The 
shaded stripe above this line also represents the region 
of the factor-of-two uncertainty. This region does not 
spread below the line, because of the nature of constraint 
(|2U)l . (Here the ten per cent uncertainty associated with 
condition (|13p and the small uncertainty in the knowl- 
edge of J are neglected.) 
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The qualitative interpretation of Fig.[6]can be phrased 
as follows: 

(i) if x C Q < 0.12, then the system almost certainly 
phase separates; 

(ii) if 0.12 < x c o < 0.18, then the system likely phase 
separates; 

(hi) if 0.18 < x c o < 0.26, then the chances of phase 
separation are about 50 per cent; 

(iv) for 0.26 < x c0 < 0.36, the chances of phase sepa- 
ration arc still significant but less then 50 per cent; 

(v) finally, for 0.36 < x c q, phase separation is increas- 
ingly unlikely. 



D. Constraint on x c q from experiments 

Recent experiments by Wakimoto et al. -32] show that 
the intensity of AF correlations integrated in the fre- 
quency range — 100 meV in L^-zSr^CuC^ at x = 0.3 is 
about 10 per cent of that in La2_ K Ba K Cu04 at x — 0.125, 
which, in turn, should be significantly smaller than the 
AF intensity in the parent compound LaCuC>4 corre- 
sponding to x = 0. Similar results but with a smaller fre- 
quency integration range, — 50meV, were also reported 
previously for YBa 2 Cu306+ x by Bourges[3lj. It thus ap- 
pears that x c o w 0.3 would satisfy condition (| 13[) defining 
the critical crossover concentration. 

However, x c q can also be significantly smaller than 0.3. 
If phase separation actually takes place in cuprates, then 
the experimental magnetic signal can come from hole- 
poor regions of the phase separated state. All we can 
conclude then is that xi « 0.3, but the maximum criti- 
cal concentration x c q for the homogeneous state can be 
significantly smaller than x^- [Here x% is the right end of 
the right metastable region defined in Section III Al and 
on Fig. EJb) .] 

We propose to further constrain x c q using the phe- 
nomenology of phase separation in LaCu04 +( 5 (Fig. [7]), 
and the fact that for a phase diagram like the one shown 
in Fig. [U the temperature-dependent critical concentra- 
tion, which we denote here as x c (T) should be close 
to the right boundary between spinodal (locally unsta- 
ble) region and metastable (locally stable) region. [Note: 
x c o = x c (0)]. As shown in Appendix [ATb). x c (T) sim- 
ply coincides with the spinodal boundary in the case of 
second order phase transitions. 

Oxygen intercalated LaCuC^+a is an atypical cuprate 
family, where macroscopic phase separation actually 
takes place, supposedly because intercalated oxygen 
atoms remain mobile between Cu02 planes above tem- 
peratures about 200K. The redistribution of intercalated 
oxygen can thus screen the uncompensated charge of elec- 
tronic inhomogeneities. A popular piece of the experi- 
mental phenomenology of phase separation in LaCuC^+a 
is shown in Fig. [Jj The horizontal axis of Fig.[7jrepresents 
the value of 6 in LaCu04+a. It is tentatively assumed to 
be related to the charge carrier concentration as x = 25. 

The energetics of intercalated oxygen can be a large 
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FIG. 7: (Color online) Experimental data points for the 
phase separation diagram of LaCu04+a from Ref. [l| with 
solid straight line superimposed on the top of it to extrap- 
olate to the right boundary of the phase separation region 
for the hypothetical case of never-frozen oxygen motion. The 
nominal charge carrier concentration is x = 28. 

part of the energy balance behind the phase separation 
in LaCu04+,5. However, it is clear from our estimate in 
section IIII C 41 that the contribution to the phase sep- 
aration energy balance from AF correlations in Cu02 
planes should also be large. When several large energy 
terms closely compete for and against phase separation 
and the overall energy gain resulting from phase sepa- 
ration is small, then it is likely that the spinodal region 
ends where one of the large competing terms suddenly 
becomes small. In the present case, we know that the 
energy of AF correlations decreases steeply in the rele- 
vant doping region. Therefore, it is natural to expect 
that the spinodal region in LaCuC^+a ends where the AF 
correlations either disappear or are drastically reduced. 

It is not entirely straightforward to extract the bound- 
ary of spinodal region from experiments. For infinitely 
slow cooling, one expects that the system starts phase 
separating once the value of X2(T) (the right boundary 
of metastable region) becomes equal to the externally 
set doping level. However, for a finite cooling rate, the 
onset of phase separation may happen at lower tempera- 
tures, because it requires activation over an energy bar- 
rier between metastable homogeneous state and phase- 
separated state. One can only be certain, that, even if 
the cooling rate is too fast (but the oxygen ions remain 
mobile), phase separation should start once the charge 
carrier concentration reaches the right spinodal bound- 
ary, which we associated with the critical concentration 
x c (T). In other words: 

(T), (21) 

where x exp (T) is the right boundary of the phase sepa- 
ration range observable experimentally at temperatures, 
where intercalated oxygen ions are still mobile. The loss 



of oxygen mobility manifests itself in the freezing of the 
experimentally observable boundaries of the phase sepa- 
ration region, i.e. these boundaries in the axes of Fig. [7] 
become vertical. One can then see that in Fig. [7] there 
are only two useful experimental points for the right end 
of the phase separation region — those through which 
the straight line is drawn. This straight line is a crude 
extrapolation aiming at finding the right boundary of the 
phase separation range for the hypothetical case of never- 
frozen oxygen ions. This line crosses the horizontal axis 
at 6 = 0.065, which, according to the assumption x = 28, 
implies x exp (0) = 0.13. Given inequality l|2"Tj) and recall- 
ing that x co = x c (0), we then obtain 

x c0 < 0.13 . (22) 

It had been reported[I(| that, in fact, x < 28. This, 
however, only strengthens the case for inequality l|22jl . 

To summarize this subsection, we have discussed the 
direct observational constraint [32| x c q < 0.3 and an in- 
direct experimental evidence (Fig. [7]) for x c q < 0.13. In 
terms of subsection IIII C 41 x c q = 0.3 implies that the 
chances of phase separation are significant but less than 
50 per cent, while x c q < 0.13 indicates that phase sepa- 
ration is likely, i.e. its chances are significantly above 50 
per cent. 

E. Theoretical discussion of microscopic factors 

The studies of the t — J and Hubbard models for the 
values of parameters relevant to the cupratesjl, [H, [3, 
EE El El, 11, SI, H| S3, IH HI agree that these models 
are close to the threshold of phase separation but do not 
reach the consensus on whether phase separation actually 
takes place or not. However, the common expectation is 
that even if the t — J or Hubbard model solutions are 
slightly over the phase separation threshold, the in-plane 
Coulomb repulsion between charge carriers (not included 
in the models but present in real materials) would eas- 
ily suppress phase separation. Ivanov[23| has reported 
numerical results directly supporting this expectation. 

The above expectation appears to be in conflict with 
the conclusion of the preceeding subsection derived from 
the analysis of experiments. Below we further suggest 
that this expectation can be misleading on purely the- 
oretical grounds. We argue that the omission of the 
"medium-range" Coulomb repulsion between the charge 
carriers separated by 1-3 lattice periods leads to a larger 
value of the critical concentration x c q, where AF cor- 
relations vanish. This, in turn, weakens the tendency 
towards phase separation. 

We introduce the effect of the medium-range Coulomb 
interactions on the doping dependence of AF correlations 
in several steps as illustrated in Fig. [8l We start from 
the half-filled state, where the exchange energy per spin 
is about J/2, and then estimate x c q as the concentration 
of holes necessary to destroy this energy. Although F v (x) 
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exchange energy 




FIG. 8: (Color online) Illustrations for various scenarios for 
the doping evolution of the AF exchange energy in cuprates: 
a) independent charge carriers; b) correlated charge carriers in 
AF background; c) medium-range Coulomb repulsion between 
charge carriers added; d) Coulomb potential of dopant atoms 
added. See the discussion in the text. 

includes indirect effects of AF correlations in addition to 
just AF exchange energy, this function should become 
zero, when AF correlations disappear, which, in turn, 
happens simultaneously with AF exchange energy falling 
to zero. 

First we consider the approximation of "independent 
charge carriers" [line (a) in Fig. [8], according to which, 
each hole independently destroys AF correlations in the 
region around itself. The AF energy cost of a doped hole 
can be estimated[4l[ using the typical length of two lat- 
tice periods for the string of broken AF bonds, which 
are left behind by the hole added to the half-filled t-J 
model with t/J = 3. This means that 1 hole creates 
approximately 7 ferromagnetic bonds at energy cost J/2 
each, i.e. it destroys AF energy about |j. This sim- 
ple estimate should be reasonable for a broader class of 
situations beyond the string picture in t — J model. 

In the above approximation, the AF energy disappears 
at the doping concentration determined by condition: 

7 -Jx cQ = \j, (23) 

which means x c q w 0.14. 

The above estimate of the AF exchange energy per 
hole should hold for very small doping concentrations. As 
Step 2 of our analysis, we take into account the fact that, 
at larger concentrations, the holes would tend to hop onto 
the sites with AF bonds broken by a nearby hole. The AF 
exchange energy cost per hole is, therefore, reduced as the 
hole concentration increases, following line (b) sketched 
in Fig. [H This sketch presumably reflects what happens 
in the t-J and Hubbard models and explains why x c q 
associated with those models can easily increase to the 
values 0.4-0.5, which, in turn, makes the onset of phase 
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separation uncertain. 

At the next step, we introduce the medium-range 
Coulomb repulsion between charge carriers. This repul- 
sion prevents one hole from taking advantage of the AF 
bonds broken by another hole. Therefore, each hole again 
destroys more AF bonds. It is this factor, that may 
have been missing in the numerical work of Ivanov[23(, 
who found that the nearest neighbor Coulomb repulsion 
suppressed phase separation. In Ref. [HI], however, a 
variational ansatz was used, which did not suppress the 
chances of holes to occupy adjacent lattice sites, and, 
therefore, not surprisingly, the resulting energy balance 
opposed phase separation. 

One can estimate the Coulomb interaction con- 
stants between holes occupying adjacent lattice sites as 
700 meV. This estimate uses the high frequency dielec- 
tric constant e ~ 5. The high frequency value is taken, 
because the lifetime of the region of broken AF bonds can 
be estimated as 1/J, which is faster than the timescale of 
lattice vibrations responsible for low frequency Coulomb 
screening. (We assume the Debye frequency 60 meV/h.) 

One should also consider the possibility of holes form- 
ing long-living pairs or larger clusters. In this case, lat- 
tice would have time to respond, which means that it 
would additionally screen the interaction between holes, 
and thereby reduce the dielectric constant to e w 30[42|. 
However, even in the latter case, the repulsion between 
holes occupying nearest neighbor sites remains significant 
- about 100 meV. 

Once the medium-range Coulomb repulsion is taken 
into consideration, the plot of exchange energy as a func- 
tion of doping concentration should pass somewhere be- 
tween lines (a) and (b) in Fig. [51 which gives line (c). 
At lower concentrations, it is closer to the independent 
hole approximation, i.e. line (a), and at higher concen- 
trations, when charge carriers have no choice but to ap- 
proach each other, the plot should be closer to line (b). 
Thus, the effect of the medium range Coulomb repulsion 
is to make the AF crossover sharper, i.e. the AF negative 
curvature in the crossover region around x w 0.14 should 
become larger. 

The next level of approximation is to introduce the 
Coulomb potential of dopant ions — line (d) in Fig.[8j At 
low doping concentrations, this potential should localize 
holes around dopant ions and thereby reduce the radius 
of the region around the hole, where AF correlations are 
destroyed. Therefore, the low doping part of the AF en- 
ergy vs. doping concentration plot should have smaller 
slope than the one in the preceeding step. As the dop- 
ing concentration increases, the Coulomb potential wells 
around the dopant ions start overlapping stronger, and 
therefore, the result of the previous step is recovered. As 
is obvious from Fig. [SJ such a correction further increases 
the AF negative curvature in the region around x « 0.14. 
Quantitatively however, this correction depends on the 
localization radius around dopant ions and may be mi- 
nor, if this radius is larger than two lattice periods. The 
latter appears to be the case for the cuprates(42j. 



Finally, the Coulomb interaction would also play an 
important role in the energy balance determining what 
kind of inhomogeneous patterns emerge if the homoge- 
neous state becomes unstable, but this is the subject be- 
yond the scope of the present work. As mentioned in 
subsection IIII C 2"lt he above instability can lead, among 
others, to stripes 0], checkerboards^. g. [H, 51]), or the 
charge imbalance between adjacent Cu02 planes. 



IV. CONCLUSIONS 

We have demonstrated that focusing on the question of 
how the energy associated with magnetic or other type of 
correlations evolves with doping, leads to a useful insight 
into the factors controlling phase separation in the vicin- 
ity of phase transitions and crossovers. We have intro- 
duced a quantitative constraint on the negative curvature 
contribution to the free energy and used this constraint 
to evaluate the chances of electronic phase separation in 
the cuprates. 

Our analysis of the antiferromagnetic crossover in the 
cuprates has led us to conclude that these materials 
are realistically close to the phase separation threshold, 
and the chances that they gcncrically exhibit some form 
of phase separation within the superconducting doping 
range are quantitatively significant. In particular, the 
lasagna phase separation scenario with macroscopic re- 
gions of opposite charge lying on the top of each other 
in adjacent Cu02 planes appears to be quite viable and 
creates additional worry for the interpretation of many 
experiments. Cuprates also appear to be closer to the 
phase separation threshold than the expectations based 
on the Hubbard and t— J models suggest. This may be re- 
lated to the fact that the above models neglect Coulomb 
interaction between charge carriers separated by 1-3 lat- 
tice sites. 

The proximity to the phase separation threshold on 
either side, rather than the presence or the absence of 
phase separation as such, may be the distinguishing char- 
acteristics of the superconducting cuprates. The result- 
ing "softness" of the electronic liquid should then lead to 
significant fluctuations of charge density and electric field 
at long wavelengths and relatively low frequencies [45|. 

This work was supported in part by the National Sci- 
ence Foundation through DMR-0404781. 
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APPENDIX A: CLASSIFICATION OF 
HOMOGENEOUS SCENARIOS BEHIND PHASE 
SEPARATION 

a. First order phase transitions 

At a first order phase transition, the system switches 
from one minimum of free energy to another, not contin- 
uously connected to the first one in the space of macro- 
scopic variables of the system. The values of free energy 
of the minima should have different dependencies on the 
carrier concentration. The switching between the two 
minima takes place at the critical carrier concentration, 
where the two free energies coincide (see Fig. [3]). As 
obvious from Fig. [3J this guarantees that the resulting 
concentration dependence of the minimal free energy has 
an upward cusp at x = x c o- The upward cusp implies 
infinite negative second derivative at x — x c q, which, in 
turn, guarantees that the system phase separates around 
x C Q [at least on nanoscale], no matter what the value of 
Coulomb interaction is. 

Since the observation of phase separation does not im- 
mediately imply that the underlying homogeneous tran- 
sition would be of the first order, one can make a com- 
pelling case in favor of such a transition, only when 
none of the two phases can emerge from the other 
as a result of spontaneous symmetry breaking: e.g., 
for antiferromagnetic — ferromagnetic transition in the 
context of colossal magnetoresistance manganites, and 
for the low temperature orthorhombic — low temperature 
tetragonal (LTO — LTT) transition in cuprates and nick- 
elates. However, even in the above two examples, the ho- 
mogeneous transition can be split into two closely spaced 
second order transitions resulting in a "rounded" cusp in 
the concentration dependence of free energy. 



Now we consider the evolution of parameters in 
Eq. (|A2|) as a function of carrier concentration as rele- 
vant to the phase diagram shown in Fig. [TJ The ratio 
may depend on the concentration x, but near (be- 
low) the critical doping concentration x c q, the key factor 
is the dependence of T c on x, which we parameterize in 
the vicinity of x c o as 



T c = \(x c0 - x), 



(A3) 



where A is a slope parameter. Any alternative param- 
eterization of the form T c = X(x c o — xp with ( > 1 or 
< C < 1 would be more favorable to phase separation 
(see Fig. [5] and the related discussion in subsection III C[) . 
Substituting Eq. fAlij) into Eq. R2|l . we obtain 



(A4) 



F v = -K v [x c (T) 



where K n the negative curvature given by 



K„ = 



4B 



and 



x c (T) 



T 

x c0 --. 



(A5) 



(A6) 



Important for the analysis of cuprates is the property 
that, if the negative curvature (| A5|) outweighs the posi- 
tive curvature of the rest of the free energy, then one of 
the ends of the miscibility gap would coincide with x c (T). 

The value of K n given by Eq. (|A5[) can now be esti- 
mated as follows: 

The slope of the T c -vs. 
crudely approximated as 



line around x = x c q can be 



b. Second order phase transitions 

We start from the standard Landau expansion for a 
second order phase transition in terms of an order pa- 
rameter having absolute value rj (e.g. the absolute value 
of staggered magnetization). The transformational prop- 
erties of this order parameter under time reversal are such 
that, in the absence of external field, only even powers of 
r) enter the expansion for the free energy, i.e. 



F v = Arf + Bif 



(Al) 



where A and B are two expansion coefficients. The onset 
of the phase transition is determined by the sign of coef- 
ficient A, which can be parameterized as A = a(T — T c ), 
where a is a positive constant, T is temperature, and T c 
is the transition temperature. At T < T c , A is nega- 
tive and hence the minimum of free energy is reached at 



>1 



q(T e -T) 
2B 



The minimized free energy is then 



F„ 



a 2 {T c -Tf 
4B 



(A2) 



A 



T c q 

XcO ' 



(A7) 



where T c q is the critical temperature at x = 0. The value 
of a 2 /(4B) can be estimated by equating the Landau 
expansion result and the approximate microscopic value 
of the gain in the ordering energy F^x) at x = 0, T = 0: 



F n o — -r^T M 



4B' 



(A8) 



The value of F v q is assumed to be known either from 
experiments or from direct microscopic considerations. 
Thus 



Kn 



)/() 
2 ' 

1 cO 



(A9) 



Despite the crudeness of the assumptions that led to 
Eq. (|A9[) . the right-hand-side of this equation is, in fact, a 
rigorous constraint on the negative curvature for a broad 
class of situations. It is formulated in subsection lll CI and 
derived in Appendix [51 
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c. Crossovers 

Crossovers constitute, perhaps, the most general case, 
where neither of the derivatives F v {x) exhibits a jump 
(see Fig. 0|). For the purposes of this work, phase tran- 
sitions of orders higher than two can also be treated 
as crossovers. Like second order phase transitions the 
crossovers must lead to the regions of the negative curva- 
ture of F v (x) . One can make a good estimate of this cur- 
vature by approximating the crossover behavior of (x) 
by a suitable second order phase transition curve. This 
procedure is put on a firm foundation in Section Til CI and 
Appendix [B] . 



F 2 








fen 


/f(x) 
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APPENDIX B: DERIVATIONS OF THE 
CONSTRAINTS ON THE CURVATURE. 

In this Appendix, we present the rigorous solutions of 
problems labeled as B, Bl and B2, which are equivalent, 
in respective order, to problems A, Al, and A2 formu- 
lated in Section HTCl While the problems A, Al and A2 
are of physical interest to us, the alternative formulation 
allows us to avoid many unnecessary minus signs and 
factors 1/2 in the solution. Everywhere below, variables 
without subscripts, like F(x), F'(x) and F"(x), are the 
functions of the argument in the parentheses, whereas 
subscripted variables like F x , F 2 , F{, F 2 ', F' max , F£ ax , 
etc., are only the numbers and, if followed by parentheses 
should be multiplied by the expression in the parentheses. 
As in Section [II CI prime and double prime superscripts 
denote the first and the second derivatives respectively. 

Problem B: Here we consider all paths leading from x = 
xi to x = X2 > xi, such that F(x\) = Fi, F'(xi - 0) = 
and F(x2) — F 2 > Fx, and the question is: What is the 
minimum possible value of the maximum positive second 
derivative (max[F" (x)]) on a path leading from F(xi) to 
F(x 2 )7 See Fig.! 

Functions F(x) need not be analytical. Here, the in- 
finite positive second derivative F"{x) would disqualify 
the path, while infinite negative second derivative may 
be present in the optimal path. With the above reserva- 
tions, we will formally treat paths F(x) as if they were 
analytical functions. 

Solution of Problem B. 

Let us denote the maximum values of F'(x) and F"(x) 
in the interval [a^a^] as F' max and F^ ax , respectively 
and the value of x, where F'(x) reaches its maximum, as 
x m . We now express F 2 as follows 



Fn 



Fi_+ F'(x) dx+ F'(x) dx 



(Bl) 



Let us now denote the maximum value of F"(x) within 
the first interval [xi,x m ] as f^ 01 . Since F'(xi) = 0, 
the value of F'(x) within the above interval can be con- 



FIG. 9: (Color online) Illustration for the boundary condi- 
tions of problem B. 



strained as follows 

F'(x)= f F"{x') dx' <F^ ax {x- Xl ). 



(B2) 



(B3) 



(B4) 



For the second integral in Eq. ljBlj) . we use the bound 



Inequality (|B2[) has two consequences 

F m ax — Fmax( x rn ~ X\) 

and 

1 F'(x) dx<^F^ ax (x m - Xi y. 



F'(x) dx < F' max (x 2 - x m ) 



(B5) 



Now, substituting the constraint (|B3|) for F' nax into in- 
equality (|B5|) and combining the result with Eq. (|F31[) and 
inequality (|B4|) , we obtain 



F 2 <F 1 + ~F^ ax [(x 2 - x,f - (x 2 - x m f] (B6) 

The expression in the square brackets above is always 
positive in the interval [x l5 x 2 ] and has maximum value, 
when x m = x 2 . Taking this into account together with 
the fact that F^ ax > F^ ax , we, finally, obtain 



F 



> 



2(F 2 - Fx) 
(x 2 - Xi) 2 ' 



(B7) 



The right-hand side of the above inequality, actually, rep- 
resents a possible value of the maximum curvature, which 
corresponds to parabola 



F(x)=F 1+ . F2 F \ 2 (x- Xl f. 

(x 2 - XlY 



(B8) 
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Hence the right-hand side of inequality (|B7|) , is the mini- 
mum possible value of F^ ax . In terms of Problem A, this 
result implies formula (|11| . The A-problem counterpart 
of parabola (|B8[) is shown in FigO 



Now we generalize Problem B to Problem Bl (equiva- 
lent to Problem Al from Section Hi Cj) . Problem Bl has 
all conditions of Problem B and, in addition: F'{xx) = 
F{ . All other notations will be the same as in the solution 
of Problem B. 

Solution of Problem Bl 

One can reduce Problem Bl for F(x) to Problem B for 
function 



G(x) = F(x) - Fl(x - xx). 



(B9) 



In this case, G{xx) = Gx = Fx, and G(x 2 ) = G 2 = 
F 2 — F{(x 2 — Xx) The curvatures of G(x) and F(x) are, 
obviously, the same everywhere . If G 2 > Gx, or, equiv- 
alently, 



F{> 



F 2 — F\ 
x 2 — xx 



(BIO) 



then inequality (|B7[) translates into 



F > 

max — 



2{G 2 -Gx) 2[F 2 - F 1 - F{(x 2 - xx)} 



(x 2 - Xi) 2 



(x 2 - xi) 2 



(BH) 

If condition (|B10|) is not satisfied, then there are paths 
between xx and x 2 , which have no regions of positive 
curvature. Inequality (|B11[) is equivalent to Eq. (TT2)) for 
Problem Al. 

Let us now turn to Problem B2, which has condition 
F'(x 2 ) = F 2 in addition to all other conditions of Prob- 
lem Bl. 

Let us further assume that F'(xx) — 0. (If F'(xx) ^ 0, 
then one can perform transformation (|B9|) and bring the 
problem to the condition desired.) 

We first consider the situation (Case I), when 



F' 2 < 



2(F 2 -F 1 ) 
x 2 - xx 



(B12) 



The right-hand side of the above inequality is the first 
derivative at x — x 2 of the optimal parabolic path (|B8|) 
achieving the minimum of F^ ax in the absence of con- 
straint on F'(x 2 ). 

Now we make the following two remarks: 

1. The minimum of F^ ax in the present case (with 
constraint on F'(x 2 )) cannot be smaller than the min- 
imum obtained for the case without the constraint on 
F'(x 2 ). Therefore, if we construct a path, whose F^ ax 
can be made arbitrarily close to the minimum F^ ax for 
the constraint-free case, then the constraint-free mini- 
mum will simultaneously be the minimum for the present 
case. 

2. Possible paths can have points of infinite negative 
second derivative (see the earlier discussion). 



The example anticipated in the first remark, indeed 
exists. The required path consists of a parabola, which 
approaches arbitrarily close from above to the optimal 
parabola (|B8[) in the case without constraint on F'{x 2 ). 
When the former parabola crosses the straight line de- 
scribed by function F(x) = F 2 + F 2 (x — x 2 ), the path 
just switches to that straight line thus exhibiting infinite 
negative second derivative. The maximum positive sec- 
ond derivative on such a path is that of its parabolic part, 
which, therefore, can have value arbitrarily close to that 
of the optimal parabola (|B8[) . Thus Case I of the present 
problem with F'(xi) = results in condition IB7I 

Now we turn to Case II, corresponding to condition 



2{F 2 -F 1 ) 
x 2 — Xx 



(B13) 



In this case, the straightforward generalization of the so- 
lution for Case I would not work, because the switching 
between the parabola and the straight line would have to 
exhibit infinite positive second derivative, which would 
be contrary to our task of minimizing F^ ax . 
Instead we define a new function 



H{x)=F{x)-F! 2 (x-x 2 ), 



(B14) 



which everywhere has the same second derivative as 
F(x). The new function is constrained by the following 
conditions: 



H{x l ) = H 1 = F 1 +F! 2 {x 2 -xx), 



H'(xx) = H[ 



H(x 2 )=H 2 =F 2 , 



H'{x 2 ) = H! 2 = 0. 



(B15) 



(B16) 



(B17) 



(B18) 



Conceptually, the problem for H(x) is a mirror reflection 
of our original problem for F{x), and, moreover, Case II 
of the problem for F{x) corresponds to Case I for H(x). 
The latter is defined by the condition 



\H[\ 



< 



2(ffi - H 2 ) 

x 2 - Xx 



(B19) 



which is, indeed, fulfilled given inequality (|B13[) and 
Eqs. (IB15llB17j) . This means that the path minimizing 
the maximum positive curvature for both H(x) and F(x) 
is a parabola, which matches the constraint on the first 
derivative at x = x 2 and switches to a straight line near 
x = Xx- It is characterized by 



F" 

max 



2(F 1 ~ F 2 + F^{x 2 - Xl )) 
(x 2 -xx) 2 



(B20) 



We conclude this Appendix by giving explicit expres- 
sions for the results of Problem B2 for the case, when 
F'(xx) + and F\x 2 ) ^ 0: 
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1. If F'(xi) > { x _ x and either there is no constraint 

on F'(x2), or F' (2:2) < ^ F2 ~ Fl ^ , then the paths between 
X\ and X2 can have no regions of positive curvature. 
Otherwise: 

2a. If either F'( Xl ) + F'{x 2 ) < 2{ x *I Xl l} or F'{x 2 ) is 
unknown, then the minimum value of F^ ax is given by 
the right-hand side of inequality (|B1 1|) . 

2b. If F'(x 2 ) is known and F'(xi)+F'(x 2 ) > 2( X *I X ' ] , 
then the minimum of F^ ax is given by Eq. (|B20p . 



In both cases 2a and 2b, the minimal values of F" lax 
correspond to parabolas passing through points (xi,Fi) 
and (x2,F2). In case 2a, the optimal parabola matches 
the given value of F'(xi), while in case 2b it matches 
F'(x 2 ). 

In the end of Section III CI the above results are reex- 
pressed for the negative curvature in terms of Problem 
A2. 
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